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This paper presents the first direct measurement of the B meson differen- 
tial cross section, da/dpT, in pp collisions at =1.8 TeV using a sample of 
19.3 ±0.7 pb -1 accumulated by the Collider Detector at Fermilab (CDF). The 
cross section is measured in the central rapidity region \y\ < 1 for pt{B) > 6.0 
GeV/c by fully reconstructing the B meson decays B + — > J/ipK + and 
B° -» J/%l)K*° (892), where J/if> -> fi+/jr and K*° -> K+tt'. A compari- 
son is made to the theoretical QCD prediction calculated at next-to-leading 
order. 

PACS numbers: 13.25.Hw, 13.85.Ni, 14.40.Nd 

The prediction of QCD theory was in agreement with the first b quark cross section 
measurements in pp collisions which were made by UAl at \fs= 630 GeV However, 
measurements at \fs= 1.8 TeV made more recently at the Collider Detector at Fermilab 
(CDF) have cast doubt on whether QCD correctly predicts either the absolute rate 

or the shape of the transverse momentum (p?) distribution. We present the first direct 
measurement of the B meson differential cross section, da/dpT, in hadronic collisions by 
measuring the mass and momentum of the B mesons decaying into exclusive final states. 
This result is a more precise test of the QCD prediction than those measurements that rely 
on model-dependent procedures to infer the b quark content and momentum distribution 
from inclusive samples. The data sample used represents 19.3 ± 0.7 pb _1 collected by CDF 
during the 1992-93 run. B mesons are reconstructed via the decays B + — > J/ipK + and 
B° — > J/tpK*°(892), with J/ip — > fi + fi~ and K*° — > K + n~ and their charge conjugates. 

Detailed descriptions of the CDF detector have been provided elsewhere ||. The com- 
ponents relevant to this analysis are briefly described here. The z-axis of the detector 
coordinate system is along the beam direction. The Central Tracking Chamber (CTC) is a 
drift chamber in a 1.4T axial magnetic field, consisting of nine superlayers, four of which 
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give stereo information. A particle must have pseudorapidity \rj\ < 1 to pass through all 
nine superlayers, which defines the rapidity range used for the cross section measurement. 

A Silicon Vertex Detector (SVX) provides high-resolution r-0 tracking information near 
the interaction region ||. The SVX detector is 51 cm long and consists of four layers 
of silicon microstrip detectors with an innermost radius of 2.9 cm. Pattern recognition is 
done by extrapolating tracks from the CTC. The longitudinal distribution of the primary 
vertices is such that SVX information is available for about 60% of all tracks. If SVX 
information is added to a track, the momentum resolution, Spt/pt, improves from ~0.002pr 
to [(0.0009p r ) 2 + (0.0066) 2 ] 1 / 2 . Surrounding the CTC are electromagnetic and hadronic 
calorimeters, outside of which are the central muon chambers, segmented into 72 modules 
which provide about 85% coverage in azimuth in the pseudorapidity range \r]\ < 0.6. 

The selection of B candidates begins by identifying J/ijj candidates that decay to two 
muons. There are three levels of trigger requirements that must be satisfied for a muon pair 
to be included in the J / ip data sample. At the first trigger level, the muons must have been 
detected in the central muon chambers and pass a minimal transverse momentum require- 
ment of ~1.4 GeV/c. Prompt muons with lower momenta range out in the calorimeters. 
The efficiency for this trigger is 90% at px = 3.8 GeV/c. At the second trigger level, at least 
one of the muon-chamber tracks must match a track found in the CTC by a hardware track 
processor. The efficiency for the track processor rises from 10% at 2.3 GeV/c to 90% at 3.4 
GeV/c. At the third(software) trigger level, detailed CTC pattern recognition and tracking 
are done, and the dimuon invariant mass is required to be within 300 MeV/c 2 of the J/ip 
mass. 

To improve the purity of the J/ip sample, the CTC track is required to match its associ- 
ated muon chamber track to within 3a in r-0 and 3.5a in z. To match the trigger thresholds, 
each muon is required to have px > 1-8 GeV/c, and at least one muon is required to have 
Pt > 2.8 GeV/c. The fitted number of J/ip reconstructed after all the requirements have 
been imposed is about 45,000. The dimuon invariant mass is required to be within 4cr of the 
known J /if) mass @. The K + and K*° candidates are required to have pt > 1-25 GeV/c; 
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furthermore, each K*° decay product is required to have px > 400 MeV/c to ensure reliable 
tracking. All pairs of oppositely charged particle tracks are considered to be candidates for 
the K*° decay products. The Ktt invariant mass is required to be within 50 MeV/c 2 of the 
K*° mass. Unlike combinatoric background, the invariant mass reconstructed for an actual 
B° meson peaks at the B° mass when the kaon and pion mass assignments are exchanged. 
To avoid double counting, only the K*° candidate with the mass closest to the K*° mass is 
used. 

We find the S-candidate mass and momentum subject to the constraint that the decay 
tracks come from a common vertex and the invariant mass of the dimuon tracks is equal the 
J/ijj mass. We require the confidence level of the fit to be greater than 0.5%. The transverse 
momentum for each B candidate is required to be greater than 6.0 GeV/c. The proper 
decay length, cr = L xy rriB/pT, is calculated, where L xy is the projection of the B vertex 
displacement onto the B transverse momentum. About 75% of the background that occurs 
when a prompt J /if) is combined with other tracks from the primary vertex is removed by 
requiring cr to be greater than 100 \im. 

The B candidates are divided into subsamples in p T ranges 6-9, 9-12, 12-15, and > 15 
GeV/c. Isospin symmetry in fragmentation, i.e. equal production of B + and B°, is expected 
to hold for central production |§. To make the best measurement of da/dpT, the B^ and 
B° data samples are combined. The resulting average B meson cross section as well as 
the separate B + cross section are shown in Table I. The remainder of the paper describes 
the analysis for the average cross section in detail. For each p? range, the B + and B° 
invariant mass distributions are fit simultaneously, using an unbinned likelihood method. 
The relative numbers of signal events for each pt bin are constrained to be proportional to 
the relative reconstruction efficiencies, e(B°)/e(B + ), which are shown in Table |. The mass 
is set to 5.274 GeV/c, which is the value found from a fit to the entire data sample. The 
mass range below 5.2 GeV/c 2 is excluded from the fits since it can include contributions 
from higher multiplicity B decay modes. The slope of the background is constrained to that 
found from a sample of candidate events that fail the vertex \ 2 confidence level requirement. 
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The B + and B° invariant mass distributions for each momentum range are shown in Fig. [I], 
and the fitted numbers of events and the statistical uncertainties are given in Table I. The 
systematic uncertainty in the fitting method is determined to be 9.2% by taking into account 
the uncertainties in the mass, background slopes, and relative reconstruction efficiencies. 
The S-meson differential cross section is calculated from the equation, 

MB) _ N 
dp T 2- C- A-e- F ■ Ap T { ' 

where N is the number of events observed, L is the integrated luminosity, A is the detector 
acceptance (including the efficiency of the kinematic cuts), e is the combined tracking and 
track-matching efficiency, F is the branching fraction, and Apt is the width of the pt bin. 
The factor of 1/2 is included because decays involving both B and B mesons have been 
reconstructed, but the quoted cross sections are for B mesons only. 

A Monte Carlo simulation employing the next-to-leading order QCD calculation || with 
renormalization scale fi = \/m b 2 + p T 2 , where the 6-quark mass, m&, is equal to 4.75 GeV/c 2 ; 
the MRS-Do proton structure functions flQ" |; the Peterson parameterization |Tl[ for fragmen- 



tation, using a value of the fragmentation parameter of 0.006; and a detector and J ftp trigger 
simulation was used to determine the acceptance, shown in Table I. 

Online (third trigger level) and offline tracking efficiencies are 97 ± 2% and 98.5 ± 1.4%, 
respectively. The efficiency of the matching requirement between the CTC track and the 
muon chamber track segment is 98.7±1%. Product branching fractions of (6.55±1.01) x 10 -5 



and (6.67 ± 1.21) x 10~ 5 pf were used for the B + -> J/ipK + and B° -> J/^K*° decays, 



which include the J/ip —> fi + fi~ and K*° — > K + it~ branching ratios. These branching ratios 
are included in the relative reconstruction efficiencies given in Table [p. 

Varying the 6-quark mass and the QCD renormalization scale used in the Monte Carlo 
simulation within their uncertainties changes the calculated acceptance by ±2% The 
systematic uncertainty in the J/ip efficiency due to the trigger parameterization was de- 
termined to be 5.5%. Additionally, a systematic uncertainty of 4% is associated with the 
reconstruction of kaons that decay inside the CTC volume. The efficiency of the 100 fim 
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cut on ct has an uncertainty of 4%, due to its dependence on the lifetime of the meson 
and on the ct resolution ||13| , which varies from 50 to 300 fim, depending on whether or not 
SVX information was available for the vertex fit. The \ 2 requirement on the common- vertex 
constraint has an additional 1% uncertainty, which was determined from the fraction of J/ip 
that fail such a requirement. Varying the polarization of the B° — > J/ipK*° decay products 
in the range (80±10)% [[Uj] changes the calculated acceptance by ±5.7%. Combining these 
values in quadrature, the reconstruction efficiency has overall systematic uncertainties of 
11.1% for the B + decay and 13.5% for the B° decay. 

The cross section results are listed in Table | and plotted in Fig. [|, where the common 
branching ratio uncertainty of 11.9% (included in Table I) is shown separately. The differ- 
ential cross section measurement at 20 GeV/c is obtained by dividing the integrated cross 
section above 15 GeV/c by an effective bin width determined from the data to be 12.7 ± 1.3 
GeV/c. The point is plotted at the mean px value of the data points. The integrated cross- 
section above 6 GeV/c is 2.39 ± 0.32 ± 0.44 fib. The solid curve in Fig. |2| shows the B 
meson differential cross section predicted by the QCD-based Monte Carlo, while the dashed 
curves indicate the variation associated with uncertainty in the b quark mass, the fragmen- 



tation parameter, and the renormalization scale W5\. The curves include the generally used 
assumption that 75% of b quarks fragment in equal amounts to B + and B° mesons P,[i~6|. 
The visual comparison between data and theory is aided by plotting (data-QCD)/QCD on 
a linear scale, as shown in the inset to Fig. |2|. To determine the level of agreement between 
the data and the theoretical prediction, the predicted cross section is fitted to the measure- 
ments, holding the shape constant and varying the magnitude. The fit yields an overall 
scale factor of 1.9 ± 0.2 ± 0.2, with a confidence level of 20%. In conclusion, we find that 
the shape of the B meson differential cross section presented here is adequately described 
by next-to-leading order QCD, while the absolute rate is at the limits of that predicted by 
typical variations in the theoretical parameters. It will be interesting to compare these data 
to higher order calculations as they become available. 
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TABLE I. B meson differential cross section, d<r(|y| < 1.0) /dpr (nb/GeV/c). The uncertainty 
listed for the relative reconstruction efficiency, e(B°)/e(B + ), includes the branching ratio uncer- 
tainty. The uncertainties shown for the number of events are statistical only. 

(pr) Acceptance e(B )/ No. of B + Total No. Averaged 

GeV/c for B + % e {B + ) B ± Events Cross section of Events Cross section 

7.4 1.28 ± 0.01 0.52 ±0.12 51 ± 10 644 ± 126 ± 139 72 ± 12 596 ± 103 ± 106 

10.4 3.71 ± 0.04 0.58 ±0.14 31 ± 7 134 ± 30 ± 29 42 ± 9 116 ±24 ±21 

13.4 6.45 ± 0.06 0.61 ±0.15 20 ± 5 50 ±13 ±11 35 ± 7 55 ± 10 ± 10 

20.0 9.55 ± 0.10 0.71 ±0.17 24 ± 4 9.5 ±1.6 ±2.6 31 ± 6 7.2 ±1.4 ±1.8 
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FIG. 1. and B° meson invariant mass distributions for the momentum ranges (a,e) 6-9 

GeV/c, (b,f) 9-12 GeV/c, (c,g) 12-15 GeV/c, and (d,h) >15 GeV/c. 
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FIG. 2. Average B meson differential cross section compared to the QCD prediction. The 
branching ratio uncertainty, which contributes 11.9% to the cross section uncertainty, is shown 
separately. The dashed curves indicate the uncertainty in the theoretical prediction. The inset 
plot shows (data-QCD) / QCD on a linear scale. 
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